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I , Abstract. Investigations of spin squeezing in ensembles of quantum particles have 

^ ' been limited primarily to a subspace of spin fluctuations and a single spatial mode 

^^ . in high-spin and spatially extended ensembles. Here, we show that a wider range 

i(»^ I of spin-squeezing is attainable in ensembles of high-spin atoms, characterized by 

f^ ' sub-quantum-limited fluctuations in several independent planes of spin-fluctuation 

ly-v , observables. Further, considering the quantum dynamics of an / = 1 ferromagnetic 

f^ ' spinor Bose-Einstein condensate, we demonstrate theoretically that a high degree of 

^^ ■ spin squeezing is attained in multiple spatial modes of a spatially extended quantum 

. . , field, and that such squeezing can be extracted from spatially resolved measurements 

^ ^ ' of magnetization and nematicity, i.e. the vector and quadrupole magnetic moments, 

^^ . of the quantum gas. Taking into account several experimental limitations, we predict 

JH I that the variance of the atomic magnetization and nematicity may be reduced as far 

■ - - ' as 20 dB below the standard quantum limits. 
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Spin squeezing describes a form of entanglement in a many-spin system for which 
the variance in measurements of particular collective spin variables is smaller than the 
minimum variance achievable in uncorrelated systems [H E]. Aside from representing 
a theoretically tractable form of entanglement in many-body quantum systems, spin 
squeezing also promises to improve measurement precision in a variety of applications 
such as chronometry, magnetometry, and atom interferometry. The recent achievement 
of metrologically beneficial pseudo-spin squeezing in cold atomic gases, on the Zeeman- 
insensitive hyperfine transitions utilized in atomic clocks, is among the first realizations 
of this promise [31 H] . 

The conventional form of spin squeezing involves the collective vector spin operator, 
F, which is the equal-weight sum of the vector spin operators for all of the ensemble's 
constituent particles. Given an average ensemble spin oriented along the z axis, the 
measurement variance of orthogonal spin projections obey an uncertainty relation, 
{{AF^y) {{AFyY) > (F;j)^/4. Ensemble states of uncorrelated particles are limited 
by the standard quantum limited variance, {{AF±Y)sql = I (-^2)1/2, achieved by 
partitioning the minimum measurement uncertainties equally between the two vector 
spin projections. 

This form of vector operator spin squeezing describes present experiments not only 
on pseudo-spin 1/2 atomic gases, i.e. where each atom is restricted to just two internal 
[3], H] or external [5] states, but also on gases of higher-spin atoms (e.g. the / = 4 
hyperfine-state manifold of Cs [6]). However, the description of higher-spin atoms by 
their vector spin alone is incomplete. On the one hand, light-atom interactions that 
affect higher-order spin moments may impair efforts to create vector-spin squeezing 
in high-spin atomic gases [71 Ej. On the other hand, the deliberate addressing and 
measurement of higher-order spin moments [9] offers benefits to certain metrological 
applications, e.g. eliminating directional dependencies of atomic magnetometers pQ] . 

Here, we consider theoretically the characterization and preparation of spin- 
squeezed states in / = 1 spinor Bose gases. We go beyond existing treatments of 
spin squeezing by evaluating such squeezing in an ensemble that is explicitly higher- 
dimensional not only in its spin degrees of freedom, but also in its spatial degrees of 
freedom. Our investigation is motivated by two experimental results: the controlled 
amplification of spin noise in ^ '^Rb quantum fluids [HI [121 [13] , and the use of a spinor 
Bose- Einstein condensate as a spatially resolving magnetometer [H]. As we discuss 
herein, these two results, respectively, provide a complete means for inducing squeezing 
of both the spin vector and quadrupole moments of the / = 1 quantum gas, and also a 
means of applying the subsequent spatially extended non-classical quantum fluid toward 
spatially resolved magnetometry with sub shot-noise sensitivity. 

1. Single-mode squeezing for arbitrary spin 

We begin by introducing a generahzed treatment of spin squeezing that applies equally 
to ensembles of spin-1/2 and also of higher-spin particles. Such squeezing will be defined 
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with respect to a coherent spin state [2], in which all particles are prepared in an identical 
single-particle spin-/ state \ip). Consider an orthogonal basis of the single-particle spin 
space that includes the state {ip) and also the 2/ states |^a). We now define many- 
body observables MJ,'^ = E„ i\Q{^\ + \i^){L\)o. and mP = E^HlDii^l - \^){L\)^ 
where the summation index a is taken over all particles in the iV-particle ensemble. 
These spin fiuctuation operators describe the subspace of spin excitations orthogonal 
to the coherent spin state. With respect to this uncorrected state, these observables 
satisfy the commutation relations {[Ma , M^ ]) = 2iN6a,b- Thus, measurements of each 
pair Ma"'^' independently obey an uncertainty relation {{AMa )'^) {(AMaY) > N'^. 
Following the treatment of squeezing for vector spin operators, we identify one mode of 
spin fiuctuations, labeled by a, to be spin squeezed when the measurement variance for 
one quadrature of the Ma -Ma plane is below the standard quantum limit, A^. 

Applying this treatment to an ensemble of spin-1/2 particles, we recall that any 
coherent spin state can be written as the m = +1/2 eigenstate of a projection of the 
dimensionless vector spin. Taking that projection to lie along z, we identify the one 
mode of spin fiuctuation operators as the Pauli operators a^ and ay. More generally, 
for an ensemble of spin-/ particles prepared in the {ip) = {mz = /) eigenstate of the Fz 
spin operator, the spin fiuctuation operators defined by the state |^i) = \niz = / — 1) 
act on I'ip) as a/2// times the spin vector operators F^ and Fy. For this one mode of 
spin-fiuctuation operators, the commutation relation matches that of the spin-vector 
operators. Thus, we recover the conventional description of vector-spin squeezing as 
pertaining to just one mode of spin fiuctuations atop a maximum-spin coherent spin 
state. 

To illustrate further the possibility of squeezing in several spin fiuctuation modes, we 
consider spin squeezing of spin-1 particles prepared initially in the jm^ = 0) state. Our 
treatment is facilitated by working in a polar state basis, where |0e) is the zero-eigenvalue 
state of F ■ e [15]. The pair of fiuctuation operators defined according to the state l^^^) 
{\4>y)) are identified as M^^^ = —N^z {—Nyz) and M^^^ = -\-Fy {—F^), where the former is 
a component of the quadrupole moment tensor [16j. The coherent spin state \mz = 0) is 
thus regarded as the uncorrelated vacuum state for independent spin fiuctuations in the 
Fx-Nyz and the Fy-N^z planes, and correlated states may be spin squeezed with respect 
to vector-spin or quadrupole-spin components, or linear combinations thereof. Finally, 
we note that this unconventional form of spin squeezing, about the polar \mz = 0) state, 
may be mapped directly onto the more conventional vector-spin squeezing discussed 
above by a unitary transformation that rotates the polar basis states onto the proper 
Fz eigenstates; experimentally, such generic unitary transformations upon alkali-atom 
hyperfine spin states are produced by rf-pulses and quadratic Zeeman energy shifts |17j . 

2. Spin squeezing in a spatially extended quantum field 

We now consider such multi-spin-mode squeezing in a spatially extended, ferromagnetic 
/ = 1 gaseous Bose-Einstein condensate where, as we show below, such squeezing 
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is produced naturally by spin- dependent interactions [IB]. As in recent experiments 
p!T|[T2]. these condensed atoms (at rest) are prepared initially in the common spin state 
\mz = 0). Within the volume occupied by the condensate, we define position-space 
spin-fiuctuation measurement densities 

M«(r)= (0t(r)0,(r) + 0t(r)0„(r)) (1) 

Mf (r) = * (0l(r)0.(r) - 0i(r)0,(r)) (2) 

where 0e(i") is now the position-space Bose field operator for the polar state \(pe), and 
the index a G {x, y} runs over the transverse spin polarizations. These observables 
represent components of the local nematicity (rank-2 moments of the atomic spin) 
and magnetization (rank-1 moments of the atomic spin) of the quantum gas, and 



obey the commutation relation 



r(l)^„^ /i^(2), 



2in{r)6'^{r — r')6a,b where 



Mr(r),M,^^^(r') 

ra(r) = {(f)l{r)(f)z{r)), and expectation values are taken with respect to the polar 
condensate. 

These measurement density operators are used to construct operators corresponding 
to spatially resolved measurements of particular spin fiuctuations of the quantum field. 
Defining the real measurement mode functions {^(r), B{r), C(r), D{r)}, we may define 
mode measurement operators as follows: 

<B,c.D \ _ f. ( A{v) B{v) \ ( M(i)(r) \ 
<kc,z. )~J [ ^(r) D{r) ) [ M(^)(r) ) ^'^ 

where the common polarization index is omitted for clarity. In this definition for 
^ABCD (and similarly for M^^^^), the local measurement weight is given as 
y/A'^ + B"^ and the locally measured quadrature is defined by the angle taxi~^ {B / A) . 
The commutation relation 

[M%,c,D. <kc.D]) = 2^ I dh {AD - BC) n{v) (4) 

is non zero to the extent that the two measurement operators probe non-parallel spin- 
fiuctuation quadratures on a common population of atoms. 

We see that our present treatment subsumes the commonly discussed single-mode 
description of spin squeezing in atomic ensembles. That is, the equal-weight sum of 
atomic spin operators is measured by spatial-mode operators defined with the uniform 
functions A = D = 1 and i? = C = 0. However, our treatment easily allows for 
a multi-mode description of spin squeezing and collective spin measurements of a 
spatially extended quantum field. For example, spin squeezing in a spatially resolved 
measurement of the atomic spin, e.g. in the spinor-gas magnetometer of Ref. [H], may be 
treated by defining measurement mode functions corresponding to every resolved pixel in 
the magnetization-sensitive image. Our treatment applies also to optical measurements 
sensitive to spatially varying linear combinations of magnetization and nematicity, which 
are achieved by measuring different components of the linear optical susceptibility tensor 

m. 
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3. Generation of spin squeezing in a ferromagnetic / = 1 Bose-Einstein 
condensate 

Now, let us consider the evolution of the paramagnetic condensate under the quantum- 
coherent spin dynamics produced by the spin-dependent contact interactions between 
atoms [201 El] • As in previous experiments, we consider the condensate to be prepared 
in the uncorrelated polar state and then allowed to evolve freely at a fixed value of the 
quadratic Zeeman shift q, which gives the difference between the average energy of the 
\mz = ±1) Zeeman sublevels and that of the jm^ = 0) state. Within the Bogoliubov 
approximation, the spin-dependent Hamiltonian may be written as follows: 



^spin = E / ^'r <^1 (^0 + ^i-t) <Pa + ^{^l' + <Pi) 



(5) 



where T-Lq = —h'^V^/2ni + q with m being the atomic mass. The spin-mixing interactions 
are described through T-Li^t = C2n(r) where n(r) > is the density of the condensate and 
C2 = 4:7ih'^{a2 — ao)/3m relates to the s-wave scattering lengths a/ for binary collisions 
among atoms with total spin / G {0, 2} [221 123] • For a ferromagnetic Bose gas, C2 < 0. 
Here, we restrict the spin excitations to the volume of the condensate. 

We find it convenient to consider separately the real and imaginary parts of (pa, 
i.e. we consider the fields Za{r) = 0^ + 0^ and Pai^^) = i{(pl, — <Pa), which commute as 
[Za{r), Pa{r')] = — 2i 6^{r — r'). Considering the linearized equations of motion for these 
operators, obtained from "Hspin given above, one identifies normal dynamical modes by 
defining a set of functions Cn(r), spanning the condensate volume, as solutions to the 
following differential equation: 

(Ho + 2Hint) Cn(r) = i?^Ho"'Cn(r). (6) 

It follows that E'^ and Cn(i") ^'^'^ real and that we may normalize the mode functions so 
that J d^v QT-LQ^Cm = 5n,m/2. With the aid of these functions, we define the magnon- 
excitation mode operators (omitting the polarization index) 

Z^= JdhUr)n,'Z{r) (7) 

dhUr)Piv). (8) 



These operators are found to be canonically conjugate, satisfying the equations of 
motion hZn = Pn and hPn = —E^Zn and commutation relations [Z„, Pm] = —i6n,m- 
The canonical evolution of magnon excitation mode n is akin to that of a harmonic 
oscillator with frequency a;„ = ^jE^jfi, with E'^> ^ defining stable magnon excitations 
and i?^ < defining unstable modes. 

The action of this canonical evolution on the initial fiuctuations in each normal 
mode may be discerned by examining the scaled covariance matrix 




Cn(t) = Re ~"\ Z„(t) P„(t) (9) 
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where the mode operators are scaled by their initial rms fluctuations. It is important to 
note that this scaling circularizes the initial fluctuations of the normal mode operators, 
but that these initial fluctuations do not have the minimum uncertainty. One flnds, 
however, that the fluctuations do reach their minimum product for excitations of a 
homogeneous condensate and approximate the minimum product for inhomogeneous 
condensates with densities varying only on long length scales. 

The Hermitian covariance matrix may be expressed as C„(t) = 'Rn{t)Sn{t)7l\^{t) 
with Suit) being diagonal with eigenvalues 



S^ = 1 + 2k sin^ ujnt ± 2 W K sin^ a;„t (l + k sin^ ujnt) (10) 

with K = (x + X"^-2)/4 and x = ^n(^n(0))/(^n(0))- We note that Ksiii^uJ is positive 
for both stable and unstable magnon modes. 

For stable magnon modes, the initial spin fluctuations rotate in the Zn-Pn plane 
and also experience periodic, bounded dilatations and contractions. In contrast, for 
unstable magnon modes, the spin fluctuations become progressively squeezed along one 
quadrature axis and amplifled along the other, with the squeezed/amplifled axes deflned 
by the rotation matrix TZn{t). The fluctuations and dynamical evolution of these mode 
operators relate directly to those of normalized spatial-mode measurement observables 
given as 

corresponding to Z„ and -P„, respectively, at t = 0. Spatially resolved spin measurements 
tailored to detecting the squeezed and amplifled quadratures of the nth normal mode are 
deflned as time-dependent linear combinations of M„ and M„ through the rotation 
matrix TZnit). 

3.1. Evolution of a homogeneous spinor condensate 

For a homogeneous condensate, the evolution following a quench of the initially 
paramagnetic quantum gas is expressed simply. Translational symmetry allows one to 
identify normal modes a priori as possessing Fourier components only with magnitude 
k. The normal-mode equation ([6]) is satisfled by the functions Ck,e = \/2/"K cos(k-r + 6'), 
with V being the condensate volume, and with the magnon-excitation spectrum given 
as El = (efc + q){€k + q — go) where e^ = h^k'^ /2m [20]. Unstable modes exist for q below 
a critical value go = — C2n, specifically within the range \e + q — go/2| < go/2. 

The initial fluctuations of the normal mode operators deflned with respect to 
the mode functions {Ck,6>} are given as {Zl^g) = (e^ + g)~^ and {Pig) = 1, giving 
K = q^/lAEl). The normalized covariance matrix for modes with wave vector k is given 
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Figure 1. Factor 5^"^ by which spin fluctuations are reduced below the standard 
quantum limits, plotted vs. evolution time t measured in units of h/qQ. Solid curves 
are shown for momentum-space spin-fluctuation modes with either maximum gain, 
obtained with ek + q = Qo/S (black), or lower gain, with e^ + q = 7qo/S (gray). For 
the latter case, the long-time limit for 5^"' is shown as a dashed line. Insets show the 
corresponding spin fluctuations in the Af'-^'-M^^-' plane (the Q representation) for the 
settings and times indicated by arrows. The long-time limits for the orientations of 
the squeezing and amplification quadratures are indicated in the insets by dashed gray 
lines. 
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cfc-l-q-go 

We see that the dynamically unstable momentum-space modes are progressively 
squeezed over time (figured]). The maximally unstable mode, with o;^ = —{qo/2HY, 
is accessed for q < go/2- For this mode, the orientation of the squeezing axis is fixed 
at an angle i!) = — 7r/4 in the M^^^-M^'^^ plane, and the greatest squeezing is achieved. 
Away from this maximally squeezed mode, the orientation of the squeezing/amplification 

■sinh^(|wfc|t) > 1, 



il 



axes varies in time. At long times, defined by Ksin^((X'fct) — ^^^i .^ 

spin fiuctuations at dynamically unstable wave vectors are squeezed by a factor that 

decreases exponentially in time and is given as 

^h^\ujk\^ 



5(-) 



,-2|a;fc|t 



fo 



(14) 



The orientation of the squeezed quadrature axis is given by the relation tan ■(9 
—h\ujk\/{ek + q), with ■(? subsequently varying between and — 7r/2. 
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This momentum-space treatment describes how the dynamical instabihties of a 
ferromagnetic / = 1 spinor Bose-Einstein condensate may be used as a mode-by- 
mode parametric amphfier, both for the amphfication of small spin fluctuations atop 
the polar condensate, as explored in Refs. [121 1131 121], and also for the generation 
of a spin-squeezed quantum field. This field may be manipulated further to suit a 
particular metrological application. For example, spin squeezing may first be prepared 
at a particular wave vector k by allowing the system to evolve for a given time at q = qi 
for which ul < 0. Thereafter, the quadratic shift may be adjusted to q = q2 at which 
the spin modes at wave vector k are dynamically stable. Under such stable evolution, 
the spin fluctuations undergo rotation in the M^-^^-M^"^^ plane, so that squeezing can be 
generated along any desired quadrature axis. More general control of the /c-dependent 
rotation angle could be realized by applying a fc-dependent optically induced quadratic 
Zeeman shift, e.g. by taking advantage of the band structure of an appropriately tuned 
and polarized optical lattice. 

4. Locality in spin squeezing and measurement 

It is illuminating to consider the evolution of a spatially extended / = 1 quantum 
field in position space. The spin-dependent contact interactions to which we have 
ascribed the generation of a correlated quantum field act locally; thus, we would expect 
spatial correlations to propagate at a finite rate across the spinor condensate. Yet, the 
above description for the homogeneous system (section 13.11) dealt with measurement 
observables that incorporate correlations between particles at infinite range. In contrast, 
given the local origin of spin correlations, we expect that spin squeezing should be 
effectively observed by local measurements of the collective atomic spin on a finite 
number of atoms. 

One approach to tracing the spatial evolution of spin correlations is to consider 
correlations among the position-space spin-fluctuation measurement densities, via the 
covariance matrix defined as 



C(r,r';t) = Re^ ^;;^J,W j^( AM«(r') AMP(r') )^^ . (15) 

For the uncorrelated condensate prepared at time t = 0, we have C(r, r'; 0) = Xn 5^(r — 
r'), i.e. correlations are purely local. At later times, for the initially homogeneous 
condensate, 

C(r,r';t) = -^^ J d'ke^^<^-^'^C,it). (16) 

Now, let us inquire as to the spatial range of spin-squeezing produced by spin- 
dependent contact interactions by asking whether finite- volume spatial-mode operators 
may be defined that will capture the spin squeezing produced in a homogeneous 
spinor condensate. We consider two possibilities. First, we define "quadrature-tuned" 
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measurement mode functions through their Fourier transform as 

^,^)^n^ni,> (1.) 

where f{k) is real. Now, following the evolution time t, evaluating the covariance matrix 
with respect to these spin-fluctuation observables, we obtain simply 



C,t{t) = ^^ j d'\.f\k)Su(t). 



We consider also a second, simpler, "uniform quadrature" mode function with Fourier 
transform 

where /Cmax is the wave vector for which maximal squeezing is achieved. 

Now we are in a position to consider forms of the measurement operators 
appropriate for capturing squeezing achieved by the spinor Bose-Einstein condensate 
amplifier. We examine the long-time limit (equation [T4l) . and consider amplification 
under two conditions. 

^.i. Shallow quench 

Under the condition go/2 < g < go, the maximally unstable mode occurs at zero wave 
vector. Expanding about k = 0, we approximate 



Vg(9o - q) f h\ ( {qo-2q) j ,^2 _ , . n l2 



h \mj \A^q{qo-q) ' 

We note that D, proportional to h/ni, has the units of a diffusion constant. Let us now 
define the measurement operator through 

fik)^'-^e-"^' (21) 



SO that the quantity {AD — BC) > is a 3D Gaussian with peak value of unity at the 
origin, and with a volume of f = (27rcr^)'^/^. 

Evaluating the correlation matrix, we now find fluctuations of the Mj^ ^^ j^ operator 
are squeezed with respect to their variance at t = 0. An estimate for the squeezing 
achieved using the quadrature-tuned mode function may be obtained analytically. If we 
assume that a^ > Dt, then the reduction in the variance is given approximately as 

4-) - ^^^^ (1 - ^) "'' e— .. (22) 

That is, the spin correlations produced in the quenched spinor Bose-Einstein condensate 
are indeed local, contained within a volume that grows diffusively in time. 

These approximations are confirmed by numerical calculations, results of which 
are shown in figure [2l The quadrature-tuned spin-fluctuation measurement is indeed 
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Figure 2. Degree of squeezing, S^~\ obtained by local measurements after preparing 
an uncorrelated \m2 — 0) condensate and quenching it to a regime of dynamical 
instabilities with q/qo = 3/4. S^^^ is calculated for variable a for quadrature-tuned 
measurement modes (solid lines) at evolution times tqo/h = {1,2,3,4,5}, and for 
uniform-quadrature (dashed lines) at times tqo/h — {1,3,5}. The quadrature-tuned 
measurement mode is found to capture the maximal squeezing produced in the system 
within a measurement volume, parameterized by cr, that grows only slowly with time. 
Open circles highlight values of a at which 5'*-~' is within 5% of the single-mode k — Q 
value. As shown in the inset, the dimension a at which this condition is achieved 
(open circles) is approximated at late times t according to equation B^ (solid line). 
In contrast, uniform-quadrature measurements must be made on larger volumes to 
observe the maximum attainable squeezing. 



effective in capturing the spin squeezing produced in this system within a measurement 
volume that grows only slowly with evolution time. In contrast, a uniform-quadrature 
measurement requires larger volumes to capture this squeezing. This requirement is 
explained by the finite differences between the quadrature axes along which squeezing 
is achieved at different wave vectors. The higher the degree of squeezing one wishes 
to detect, the more precisely must the measurement quadrature axes be tuned to the 
spin-squeezed axis by restriction to a narrow range of wave vectors. 



4-2. Deep quench 

For a deep quench, with q < go/2 the maximum instability occurs on a sphere in k- 
space, with radius fcmax defined by e^ -f g = go/2. About this shell, the gain is given 
approximately as 



Wfc 



2h 



h qo-2q 



m 



Qo 



I hj hjYi 



2h 



D{k~k^ 



(23) 
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Figure 3. Degree of squeezing, S'^~', obtained by local measurements following a 
deep quench to q = — 2go- Quadrature-tuned measurements taken at evolution times 
tqo/fi = {1, 2, 3, 4, 5} evince the maximal single-mode spin squeezing for mode functions 
given by equation [24] with a growing only slowly with time (solid lines) . Open circles 
highlight values of a at which S^^^ is within 5% of the single-mode k — fcmax value. As 
shown in the inset, the dimension a at which this condition is achieved (open circles) is 
approximated at late times t according to equation l25l (solid line). Uniform-quadrature 
mode functions, shown here for tq^/h — {1,3}, require much larger a to capture the 
full spin squeezing in this system. 



To capture this squeezing, we define a measurement mode through the following relation: 

/(fc) = Xe-"'(^-'='"^^)' (24) 

with X chosen so that AD — BC = 1 at r = 0. For large a, the volume defined by this 
function tends toward v = V^n^ ^/^max> ^^ ^'^^^ ^^^ effective radius of the measurement 
volume is r ~ (c/^max)"'^^^- 

Now we determine the degree of spin squeezing captured by this measurement 
operator. Repeating calculations as above, with the requirement o"^ > Dt, we have 



o(-) 

'-'qt 



-^) 



-1/2 



-qot/h 



(25) 



Again, we see that local measurements will capture the spin squeezing produced by the 
spinor Bose-Einstein condensate amplifier. Here, the correlation volume grows slower 
than diffusion, with radius scaling only as r oc t^'^. Numerical evaluations of S^~' are 
presented in figure [3l confirming the rapid convergence onto the single-mode maximal 
degree of squeezing for finite- volume, quadrature-tuned measurement modes. 
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5. Experimental considerations 

We conclude with a discussion of experimental prospects for realizing and detecting 
multi-mode spin squeezing in a quenched ferromagnetic / = 1 condensate. We focus 
on the case of ^^Rb, which has been shown experimentally |25l EGJ |27] and theoretically 
[281 [29] to be ferromagnetic, with Aa = — 1.4(3) a^ with as being the Bohr radius. 
The behaviour of condensates prepared in the jm^ = 0) state and quenched to the 
regime of dynamic instabilities by a rapid change of q have been studied using in situ 
magnetization-sensitive imaging. Pertinent to the present discussion, the amplification 
of spin fluctuations following quenches to different values of q was measured precisely 
in Ref. [12], and such amplification was found to be roughly consistent with quantum- 
limited amplification of quantum fluctuations in the initial state, with a gain, measured 
by the increase in the variance of spin fluctuations, as high as +30 dB. 

In contrast to previous studies on the parametric amplification of spin fluctuations 
[T2| [2l] , here we have considered the squeezing that may accompany such amplification. 
We assume experimental conditions similar to those of Refs. [HI [12], i.e. an optically 
trapped / = 1 ^^Rb condensate of A^ = 2.0 x 10^ atoms with peak density n = 
2.6 X 10^'^ cm~'^. At this density, the characteristic time scale for spin mixing in the spinor 
condensate is h/qo = 8 ms and the characteristic length scale is {h^ / {2mqo)Y^'^ = 1.8 /^m. 
Let us estimate the maximum observable degree of spin squeezing, taking into account 
non-linearities, atom loss, and the spatial inhomogeneity of the condensate. 

5.1. Single-mode nonlinearities and atom loss 

Thus far, we have treated spin dynamics in a quenched / = 1 spinor condensate 
exclusively through the Bogoliubov approximation (equation [5]), i.e. neglecting high- 
order terms in the spin fluctuation operators and also the depletion of the \mz = 0) 
condensate due to spin mixing. Within this linear, undepleted-pump approximation, 
we obtain the unphysical result, expressed in equation [TOl that the minimum spin- 
quadrature variance S^~^ in each unstable spin mode tends exponentially to zero. As 
in quantum-optical systems, sensible limits to squeezing are recovered by going beyond 
the linearized equations of motion and the undepleted-pump approximation. 

These limits are partly elucidated in the full dynamics of a single spin mode. We 
consider therefore the dynamics of a ferromagnetic spinor condensate in the single-mode 
approximation. Such a treatment would be appropriate for a spinor condensate spanning 
a volume v < {h'^/2mqo)^''^, realized experimentally with an atom number of nf < 1000. 
The full Hamiltonian is now given as 






C2 ,+ , ,+ , 1 C2 



^-rz-r.-ra-r.^ 2t; 



+ ^0t0,0t0, -^ (0t0^ _ 0t0^y (26) 



where the Bose operators (pe annihilate particles in polar state \4>e) and in the single 
spatial mode of the condensate. 
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In this expression, the single-polarization term oc ( 0], 0^ + 0j 0^ ) matches the 
two- axis counter-twisting Hamiltonian described in Ref. [2]. On its own, this number- 
conserving Hamiltonian produces squeezing that saturates with time due to depletion 
of the pump (the population ^J^^)? reaching the Heisenberg limit with S^'"^ oc 1/A^ 
[30] . However, the Hamiltonian contains additional terms: the externally imposed 
and pump-dependent quadratic Zeeman shifts and inter-polarization coupling. These 
nonlinear terms restrict spin squeezing to sub-Heisenberg scaling, with S*^"^ oc 1/yN . 
Nevertheless, the expected single-mode squeezing is significant. We have performed 
exact single-mode calculations for a 1000-atom sample, and find a reduction of spin 
fluctuations by as much as 20 dB. 

The single-mode limit allows us to discuss an additional limitation arising from the 
decay of atoms from the trap. Atom loss can be modeled by adding a noise source and 
a density-dependent loss rate in the coherent evolution of the Bose fields [30]. We have 
performed numerical calculations including such terms, representing quantum noise in 
the evolution of a classical field as in the truncated Wigner approximation (TWA) [31] . 
In accord with experimental observations, we assume an initial per-atom loss rate of 1 
s^^. We find squeezing in a A^ = 1000 atom sample can still be the level of -20 dB, a 
result that is not surprising given the very short spin-mixing time as compared to the 
lifetime of trapped atoms. 

5.2. Effects of inhomogeneity and multi-mode nonlinearities 

Experimental realizations of spinor Bose-Einstein condensates also differ from our 
idealized treatment in that the condensates are inhomogeneous, and, thus, the adequacy 
of a linearized normal-mode description to long-time dynamical evolution is less clear 
than in the homogeneous case. To assess the influence of mixing between normal modes 
and also atom losses for this situation, we have performed numerical simulations of 
the evolution of a quasi-one- dimensional quenched spinor Bose-Einstein condensate. 
Such a condensate is assumed sufficiently tightly confined in two radial dimensions 
that the spinor wavefunction is radially uniform, while being more weakly confined in 
the third dimension, along which the wavefunction is allowed to vary. Specifically, we 
consider a condensate with Thomas-Fermi radii of 3 and 200 /im along the radial and 
axial dimensions, respectively. The central, radially averaged condensate density Uq 
determines the spin mixing interaction strength to be 2|c2|t2o = h x 15 Hz. 

The quantum evolution of this spinor condensate is calculated using the TWA 
to simulate the effects of quantum noise and squeezing [24] . This approximation 
involves adding random noise to the initial mean-field condensate wavefunction, with 
the noise variance determined, in this case, by the one- dimensional condensate density. 
Temporal evolution is thereafter considered according to classical wave equations, which 
include the full non-linear dynamics given by spin-mixing Hamiltonian (the non-linear 
Schrodinger equation) and also the effects of atom loss (through additional noise and 
loss terms, described above). Expectation values, variances, and correlations for various 
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Figure 4. Numerical simulations of amplification and squeezing of spin fluctuations 
in a quenched, trapped spinor Bose-Einstein condensate. The evolution of quantum 
noise is simulated via the TWA. (a) The results of simulation runs, each with different 
instances of random initial noise, are shown at t = 65 ms after a quench to the regime of 
dynamical instabilities (see text for details) . One polarization of local spin fluctuation 
observables (components of magnetization (red solid line) and nematicity (blue dashed 
line)) is shown, scaled to unity for a fully spin polarized condensate at its center. Strong 
local correlation of the amplified spin fiuctuations is evident, (b) Squared values of 
the amplified (solid symbols) and squeezed (closed symbols) normal-mode quadrature 
measurement operators, evaluated for runs 1 (circles), 2 (squares) and 3 (diamonds), 
scaled to their initial statistical variance, indicate a reduction of initial spin fluctuations 
by 20 to 30 dB. 



measurement operators are quantified by repeating the calculations for many instances 
of the initial noise. These calculations were found convergent at a temporal resolution 
of 3.5 /is and spatial resolution of 0.5 /xm. 

We present results of a numerical simulation of an initially polar, trapped, ^^Rb 
/ = 1 spinor condensate quenched at t = to a quadratic Zeeman shift of g = /i x 2.5 Hz. 
The results for several numerical runs, i.e. for three random instances of initial noise as 
evaluated according to the TWA, are presented in figure HIi. Dynamical instabilities 
following the quench to low q result in amplified fluctuations in the transverse 
components of the magnetization and the related components of the nematicity; only 
one polarization of the observables, F^ and Nyz, is shown in the figure. The strong local 
correlation between these observables is clearly evident, reflecting the i) ~ ir/A angle of 
the amplified quadrature axis for the maximally unstable normal modes. 

Less evident in these position-space plots is the sharp reduction in the initial 
fluctuations of the normal-mode measurement observables. To reveal this reduction, we 
identify the normal dynamical modes based on the initial condensate distribution. The 
temporal gain determined for the flrst 41 unstable modes is shown in the inset of flgure 



Spin squeezing of high-spin, spatially extended quantum fields 



15 



mode index 






1 


10 

1 


20 

1 


30 

1 


40 

1 


f" 40- 


octxh 


'*'*'»=^ 


'^ajci, 


*^, 




c 20- 








■^^ 


^■ 












O 




10 



1 

20 
mode index 



~~r~ 
30 



10 



.(-) 



0.1 



r 0.01 



40 



J- 0.001 



Figure 5. Suppression of spin fluctuations in a quenched, trapped spinor Bose- 
Einstein condensate, evaluated for normal dynamical modes (labeled by mode index 
n) at several times of evolution (labeled on left of graph). The calculated early-time 
gain of these unstable modes is shown in the inset. The measurement variance in 
the squeezed quadrature mode S^~^ diminishes exponentially in time, saturating at an 
evolution time of i = 65 ms to a minimum value of 3.6 x lO"'^ for the highest-gain 
mode, and then growing due to saturation and nonlinear effects. 



El sorted in decreasing order of temporal gain (|a;„|) and labeled by an integer mode 
index. Normalized spatial-mode measurement operators are then formed as defined in 
equations [H] and [121 and the values of these operators at different evolution times t are 
evaluated for an ensemble of 600 numerical runs. From these values, we evaluate the 
scaled covariance matrices C„ numerically and determine therefrom the quadrature axes 
that are amplified and squeezed, and also the measurement variance along those axes. 

The squares of the measurement outcomes along the amplified and squeezed 
quadrature axes for the three simulation runs discussed above are shown in figure Hb- 
Aside from pronounced variations according to the initial random values assigned to the 
quantum noise in these simulations, these data indicate a suppression of spin fluctuations 
between 20 and 30 dB. 

The ensemble averaged reduction S^^^ in the initial spin fluctuations, determined for 
each mode n and at different evolution times, is shown in figure |5l We find the temporal 
evolution of the dynamically unstable modes causes the measurement variance of the 
squeezed quadrature operators to diminish exponentially in time for evolution times 
up to around 60 ms, after which the variance saturates and then begins growing. The 
minimum variance of around 3.6 x 10~^ or -24 dB is achieved for the unstable modes 
with highest gain (lowest mode index) at t = 65 ms. This maximum degree of squeezing 
matches well with the estimated effects of nonlinearities and atom loss presented above. 
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6. Conclusion 

We have presented a formalism for the evaluation of spin squeezing in ensembles 
of high-spin atoms. This formalism identifies several independent planes of spin 
fluctuations atop an arbitrary coherent spin state, revealing a greater resource for 
applications of spin squeezing in quantum information and metrology than for ensembles 
of spin- 1/2 particles. Moreover, we developed an understanding of spin squeezing 
in spatially extended quantum fields, where squeezing in many spatial modes can be 
separately generated and measured. This formalism is applicable to many contemporary 
experiments on spin squeezing of atomic gases or of other distributed quantum objects. 

Applying this formalism, we consider a source of spin squeezing in an extended 
spinor Bose-Einstein condensate. Referring to recent experiments studying the evolution 
of ®^Rb condensates quenched to a regime of dynamical instabilities [HI [111 [13] , we show 
that the strong amplification of spin fluctuations demonstrated in those experiments is 
accompanied also by spin squeezing, i.e. that the dynamically unstable normal modes act 
as quantum parametric amplifiers of spin fluctuations. Treating experimental realities 
such as atom losses and effects of confinement, we apply analytic and numerical tools to 
determine that measurable spin squeezing better than 20 dB may be achieved in such 
systems. 

This treatment is relevant also to the question of symmetry breaking in a quantum 
many-body system. An initially polar-state spinor Bose-Einstein condensate preserves 
the 5*0(2) rotational symmetry about the alignment axis. The quantum quench 
considered in this work, a rapid reduction in the quadratic Zeeman shift to induce 
dynamical instabilities and parametric amplification of spin fluctuations, traverses 
the atomic system across a symmetry-breaking quantum phase transition, where 
the spin-dependent contact interactions favor broken symmetry states of maximum 
magnetization [H]. In a closed noiseless quantum system, bona fide symmetry breaking 
does not occur; rather, the system evolves coherently to a symmetric superposition 
of broken-symmetry states. The spin-squeezed states described in this work indeed 
represent such a coherent superposition. 

The metrological significance of spin squeezed quantum fields, and also the question 
of symmetry breaking in quantum systems, provide strong motivation for the detection 
and characterization of squeezing in quenched spinor Bose gases. For this purpose, 
spatially resolved measurements of the magnetization of optically trapped condensates 
have been demonstrated [H], and optical means of measuring components of both 
the magnetization and the nematicity have been suggested [19]. However, sensitivity 
below the atomic shot noise level, required for the direct detection of squeezing, is 
experimentally challenging and was not achieved in those measurements. 

Alternately, we propose that one could use the condensate itself as a pre- 
measurement amplifier. Following the generation of inhomogeneous spin squeezing with 
the quadratic shift at gi, one may quickly increase the quadratic Zeeman shift. The 
previously unstable dynamical modes now evolve stably as the squeezed and amplified 
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spin-fluctuation quadratures rotate in their quadrature spaces. After tlieir rotation, 
tlie spin fluctuations may be again parametrically amplified by returning the quadratic 
shift to qi. Spin squeezing produced by the first stage of amplification would result in a 
reduction of spin fiuctuations observed after the second amplification stage, as compared 
to a single-stage amplification of vacuum spin noise. We have performed TWA-based 
numerical simulations which confirm that this technique is applicable to trapped spinor 
Bose-Einstein condensates at experimentally accessible settings. 
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